Mon. Not. R. Astron. Soc. 000, [Tfllip (2005) Printed 2 February 2008 (MN HT E X style file v2.2) 



The Las Campanas Infra-red Survey. V. Keck 
Spectroscopy of a large sample of Extremely Red Objects 

M. Doherty, 1 * A.J. Bunker, 2 R.S. Ellis 3 and P.J. McCarthy, 4 

1 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA, U.K. 
2 School of Physics, University of Exeter, Stacker Road, Exeter, EX44QL, U.K. 

3 California Institute of Technology, Astronomy, Mail Stop 105-24, Pasadena, CA 91125, U.S.A. 

4 Observatories of the Carnegie Institute of Washington, Santa Barbara Street, Pasadena, CA 91101, U.S.A. 



Released 2005 Xxxxx XX 



ABSTRACT 

We present deep Keck spectroscopy, using the DEIMOS and LRIS spectrographs, 
of a large and representative sample of 67 "Extremely Red Objects" (EROs) to H = 
20.5 in three fields (SSA22, Chandra Deep Field South and the NTT Deep Field) 
drawn from the Las Campanas Infrared Survey. Using the colour cut (J — H) > 3.0 
(Vega magnitudes) adopted in earlier papers in this series, we verify the efficiency of 
this selection for locating and studying distant old sources. Spectroscopic redshifts 
are determined for 44 sources, of which only two are contaminating low mass stars. 
When allowance is made for incompleteness, the spectroscopic redshift distribution 
closely matches that predicted earlier on the basis of photometric data. Our spectra 
are of sufficient quality that we can address the important question of the nature 
and homogeneity of the z >0.8 ERO population. A dominant old stellar population 
is inferred for 75% of our spectroscopic sample; a higher fraction that than seen in 
smaller, less-complete samples with broader photometric selection criteria (e.g. R—K). 
However, only 28% have spectra with no evidence of recent star formation activity, 
such as would be expected for a strictly passively-evolving population. More than 
^30% of our absorption line spectra are of the 'E+A' type with prominent Balmcr 
absorption consistent, on average, with mass growth of 5-15% in the past Gyr. We 
use our spectroscopic redshifts to improve earlier estimates of the spatial clustering 
of this population as well as to understand the significant field-to-field variation. Our 
spectroscopy enables us to pinpoint a filamentary structure at z = 1.22 in the Chandra 
Deep Field South. Overall, our study suggests that the bulk of the ERO population 
is an established population of clustered massive galaxies undergoing intermittent 
activity consistent with continued growth over the redshift interval 0.8< z <1.6 

Key words: galaxies: evolution - galaxies: stellar content - cosmology: large-scale 
structure of universe - galaxies: distances and redshifts 



1 INTRODUCTION 

One of the major challenges facing the hierarchical models 
of galaxy formation (e.g. Somerville et al. 2004) is posed 
by the large number of massive galaxies with apparently 
well-established stellar populations observed at high red- 
shift (Glazebrook et al. 2004, Cimatti et al. 2004). Galax- 
ies with evolved stellar populations at moderate redshifts 
(z ~ 1 — 2) can be characterized by very red optical-near- 
infrared colours. The first examples were identified shortly 
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after the introduction of near-infrared imaging arrays (El- 
ston, Rieke & Rieke, 1988; McCarthy, Persson, & West 1992; 
Cowie et al. 1990). Spectroscopy of two of the Elston et al. 
(1988) red candidates revealed them to be luminous evolved 
galaxies at z ~ 0.8 (Elston et al. 1989). More distant red 
galaxies with evolved populations were later identified from 
deep radio surveys (Dunlop et al. 1996). 

The optical-to-near-infrared colour selection appears to 
yield a mixture of evolved systems and actively star-forming 
galaxies with substantial internal reddening (see McCarthy 
2004 for a review). The prototype of the latter class (ERO 
J164502+4626.4) was identified by its extremely red R — K 
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colour by Hu & Ridgway (1994). Subsequent spectroscopy 
and sub-mm observations revealed this system to be a star- 
forming galaxy at z — 1.44 (Graham & Dey 1996; Cimatti 
et al. 1998, Dey et al. 1999). 

Earlier papers in this series (McCarthy et al 2001, Firth 
et al. 2002; Chen et al. 2002) presented source counts, an- 
gular clustering and photometric redshifts for objects se- 
lected via optical-near infrared colours using the Las Cam- 
panas Infra-red Survey (LCIRS) - one of the first to use a 
panoramic near-infrared camera. This work has been com- 
plemented by independent surveys of /("-band selected sam- 
ples (e.g. Cimatti et al. 2002 - K20; Yan et al. 2003 ; Abra- 
ham et al. 2004 - GDDS). Although the abundance of EROs 
is less than that implied if all of the present-day spheroidal 
galaxies followed a passive evolutionary track from high red- 
shift, it is much higher than that predicted by current semi- 
analytical models (Firth et al. 2002; see also Cimatti et al. 
2002). 

Strong clustering has been observed for the ERO pop- 
ulation, with amplitudes up to 10 times that seen in the 
equivalent flux-limited sample (see e.g. Daddi et al. 2000, 
Firth et al. 2002). On the basis of photometric redshifts, 
Firth et al. (2002) and McCarthy et al. (2001) derived co- 
moving correlation lengths of r ~ 6 — 10/t -1 Mpc (where 
h=Ho/100 km s -1 Mpc -1 ), comparable to that for early- 
type galaxies at low redshift. Daddi et al. (2002) infer a 
similar, although somewhat larger, clustering scale on the 
basis of Passive Luminosity Evolution (PLE) models of the 
redshift distribution of the K20 sample. 

Spectroscopic data has played a key role in furthering 
our understanding of the ERO population. Dunlop et al. 
(1996) and Spinrad et al. (1997) used low resolution spec- 
troscopy to infer an age in excess of 3.5 Gyr for the radio- 
selected red galaxy LBDS 53W091 at z - 1.5. More re- 
cent intermediate dispersion spectra for co-added samples 
of fainter colour-selected objects at 1.3 < z < 2 also point 
towards early formation redshifts for a significant fraction 
of the red population (McCarthy et al. 2004; Cimatti et al. 
2004). 

The goal of this paper is to build on the earlier LCIRS 
photometric studies of EROs, using deep Keck spectroscopy. 
We investigate the redshift distribution and spectroscopic 
properties of a representative subset of the ERO population. 
In order to verify and further exploit the earlier work, our 
spectroscopic targets were selected from three of the LCIRS 
fields according to the same I — H colour cut (McCarthy et 
al. 2001). 

With the resulting spectroscopic sample we aim to: 

i) verify that the I — H colour selection is optimal for 
selecting passively evolving galaxies at z > 1, 

ii) determine the redshift distribution of the population 
and compare it with that determined photometrically, 

iii) examine the spectroscopic nature of EROs with re- 
spect to the (perhaps idealized) view that such a colour 
criterion locates a passively-evolving population of sources 
destined to emerge as present-day spheroidal galaxies. 

It should be borne in mind that spectroscopic (or mor- 
phological) studies can present an incomplete picture of the 
evolutionary history of a given population. The age of a stel- 
lar population need not necessarily coincide with the age 
of the assembled mass. In particular, a spectral classifica- 



tion indicative of an established or aged stellar population 
need not necessarily correspond to a morphological ellipti- 
cal. However, spectral diagnostics of recent activity can be 
invaluable indicators of continuing growth. Our purpose in 
this paper is largely to explore the spectroscopic homogene- 
ity of the ERO population as well as the extent to which the 
idealized picture of a passively-evolving component can be 
made to fit the data. We will examine and compare the mor- 
phological information on a similar population in the next 
paper in this series (Doherty et al. in prep.) 

A plan of the paper follows: in Section 2 we present the 
sample selection and in Section 3 we detail the Keck spectro- 
scopic observations and their reduction. Section 4 presents 
the inferred redshift distribution taking into account incom- 
pleteness effects. We analyze our results in terms of our pri- 
mary objectives (above) in Section 5. Our conclusions are 
presented in Section 6. 

Unless otherwise stated, throughout the paper we use 
the standard "concordance" cosmology of Dm = 0.3, = 
0.7, and Ho = 70kms -1 Mpc -1 . All magnitudes are on the 
Vega system. 



2 SAMPLE SELECTION 

Our selection of Extremely Red Objects is drawn from three 
(out of a possible five) LCIRS fields at high Galactic latitude 
accessible at the Keck observatory: the Chandra Deep Field 
South (CDFS; Giacconi et al. 2001 ), the NTT Deep Field 
( Arnouts et al. 1999 ) and SS A22 (part of the CFRS redshift 
survey, Lilly et al. 1995 ). 

As discussed, we employ the same colour criterion 
(J — H > 3) for EROs as that adopted in the earlier photo- 
metric study by McCarthy et al. (2001). Figure demon- 
strates that this colour cut is designed to select passive 
objects with redshifts z > 1 and starbursts with substan- 
tial dust reddening. To facilitate a reasonable success rate 
with optical spectroscopy, we adopted a magnitude limit of 
H = 20.5. At this limit, the LCIR survey is close to 100% 
complete (see discussion in Chen et al. 2002 - §4.3 and 
Figure 5) and galaxies with / — H > 3 comprise roughly 
~ 10% of the total population, with a surface density of 
fa 0.8 arcmin -2 . 

We briefly recap the photometric data which forms the 
basis of our sample. The infrared imaging was obtained 
with the CIRSI camera (Beckett et al. 1998) on the Du 
Pont 2.5 m telescope at Las Campanas Observatory between 
1998 and 2000. CIRSI has 4 HgCdTe 1024 2 infrared (IR) ar- 
rays, spaced by 0.9 array widths, with a pixel scale of 0.2". 
By observing a 2 x 2 mosaic, a contiguous 13xl3arcmin 
(170 arcmin 2 ) field was surveyed. In the case of SSA22, only 
imaging data from half the area was reduced in time to se- 
lect EROs for LRIS mask manufacture for our October 2001 
observing run. Hence the total survey area from which we 
built catalogues of EROs is 425 arcmin 2 . The photometric 
data reduction is detailed in Firth et al. (2002) and Chen 
et al. (2002). In the reduced //-band images, the typical 5cr 
detection is H — 21.5 in a 3"-diameter aperture. The optical 
imaging came from the CFHT 12k Mosaic (V-, R- & /-band 
for SSA22), the Wide Field Camera on the 2.5m Isaac New- 
ton Telescope (for NTT), and the Mosaic-II camera on the 
CTIO 4m (V-, R- I- &z-band for CDF-S). The optical and 
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Figure 1. (a) Evolution in the predicted (/ — H) colours (Vega system) for various stellar populations, produced with Bruzual & 
Chariot's (2000) synthesis code assuming a Scalo (1986) initial mass function. The elliptical (E) spectral energy distribution (solid line) 
has an age of 3 Gyr. The E+A spectrum is equivalent but viewed 100 Myr after a burst of star formation involving 5% of the stellar 
mass. The starburst galaxy assumes an exponential star formation history with an e-folding time of 1 Gyr, viewed 10 Myr after the onset 
of star formation. Reddening of E(B — V) = 0.5 and 1.0 have been applied using Calzetti's (1997) prescription. The I — H > 3 colour 
criterion (shown) is designed to select passive objects with redshifts z > 1 and highly reddened starbursts with z > 0.7. For comparison, 
the equivalent tracks for an R — K > 5 colour cut are shown in the adjacent panel (b). The essential difference evident here is that E+A 
type galaxies would not be expected to occur in an R — K > 5 selected sample. 



infrared images were registered and distortion-corrected to 
an astrometric frame determined using the Digitized Sky 
Survey. Colour selection was performed using version 2.2.1 
of the SExtractor package (Bertin & Arnouts 1996) in two- 
image mode with a fixed aperture of 3"diameter. 

Spectroscopy was performed at Keck with both the 
Low Resolution Imaging Spectrograph (LRIS; Oke et al. 
1995) and the Deep Imaging Multi-Object Spectrograph 
(DEIMOS, Faber et al. 2003; Phillips et al. 2002). The 
5' x 7' LRIS slitmask field is much smaller than a CIRSI 
tile, whereas the DEIMOS field covers 16.5' x 5' (i.e., about 
half the area surveyed in a CIRSI tile). In selecting sources 
for multi-object spectroscopy, the astrometric centre and po- 
sition angle of the masks were optimised within the LCIRS 
fields to maximize the number of candidate EROs. We dis- 
cuss later the possible biases that this may bring to our 
analyses. The minimum slit length was set to 6", and the 
selection was determined entirely by geometric constraints. 
In the event of a slit clash higher priority was given to the 
brighter source (in I). 

The fraction of EROs within the field of view of each 
spectrograph that could be incorporated onto a mask was 
roughly 50% (29/73 in NTT, 16/31 in CDFS, 30/55 in 
SSA22). The angular distribution of EROs is shown in Fig- 
ure Part of the DEIMOS mask for CDFS was taken up 
with targets for another program involving spectroscopy 
of i'— band drop-out redshift ~ 6 galaxies (see Bunker et 
al. 2003; Stanway et al. 2004), which necessitated part of 
the mask falling outside of the LCIRS field. The spectro- 
scopically surveyed area from each LCIRS image is thus: 



« 80arcmin 2 for NTT; w 60arcmin 2 for CDFS; and w 
60arcmin 2 for SSA22, a total of about 200arcmin 2 . 



3 SPECTROSCOPY 

3.1 LRIS multi-object spectroscopy 

We obtained long-slit spectroscopy of 41 EROs on the nights 
of 2001 October 22-24 UT using LRIS (Oke et al. 1995) at 
the f/15 Cassegrain focus of the 10-m Keck i telescope. Spec- 
tra were obtained for 11 EROs in CDFS (including 6 which 
overlapped with the later DEIMOS sample, Section 13.21 . 
and for 30 EROs in SSA22 (using 2 different masks). Full 
details are given in Tabled 

The LRIS red-arm detector is a Tek 2048 2 CCD with 
24 /im pixels. The angular scale is 0'.'212 pixel -1 , and the 
CCD was read out in two-amplifier mode. Observations were 
obtained using the 600 line mm -1 grating in first order 
blazed at 7500 A, producing a dispersion of 1.24 A pixel -1 . 
The reference arc lamps and sky-lines have a full width at 
half maximum (FWHM) ps 5-8 A (the spectral focus being 
best at the central wavelength). For objects that fill a 1"- 
wide slit, the velocity width of a spectrally unresolved line 
is FWHM ss 220 - 300 km s" 1 . A wider 1.5" slit was used 
for CDFS and SSA22#1, producing slightly lower spectral 
resolution (~ 400km s -1 FWHM). The grating was tilted 
to place a central wavelength of A = 8200 A on the detector 
(for slits in the centre of the mask) , sampling the wavelength 
range 7000 - 9000 A in all targets, with up to ta 500 A ei- 
ther end depending on the position of the slit on the mask. 
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Field 


Field Centre (J2000) 


Instr. 


date/time 


total exptimc 


PA (deg) 


slit width (arcsec) 


No. targets 


SSA22#1 


22:17:41.9 00:13:53 


LRIS 


22-24 Oct 2001 


16ks 


10 


1.5 


16 


SSA22#2 


22:17:45.0 00:18:45 


LRIS 


22-24 Oct 2001 


26ks 


-20 


1.0 


12 


CDFS 


03:32:26.1 -27:43:24 


LRIS 


22-24 Oct 2001 


18ks 


-25 


1.5 


ll 1 


CDFS 


03:32:35.2 -27:47:52 


DEIMOS 


8-9 Jan 2003 


20 ks 


6 


1.0 


15 2 


NTT 


12:04:10.2 -07:26:06 


DEIMOS 


8-9 Jan 2003 


19.2 ks 


-50 


1.0 


23 



Table 1. Summary of Spectroscopic Observations 



We simultaneously obtained exposures with the blue arm of 
LRIS (read in one-amplifier mode, with the 300 lines mm - 
grating blazed at 5000 A), using a dichroic beam-splitter at 
6800 A. However, due to the red nature of our objects, there 
was little or no flux at A < 6000 A. 

The total integration times are given in Table These 
were broken into individual exposures of duration 2000 s to 
enable more effective cosmic ray rejection. The telescope was 
dithered by ±3" along the slit between integrations to facil- 
itate the removal of fringing in the red and the elimination 
of bad pixels. The observations spanned an airmass range 
of 1.00 - 1.5. The seeing had a FWHM = O'.'6-l'.'O over the 
course of the night. Spectrophotometric standard stars HZ 4, 
G191B2B & Feige 110 (Massey et al. 1988; Massey & Gron- 
wall 1990) were observed at similar airmass to determine the 
sensitivity function for flux calibration. 

Each frame was bias subtracted according to values de- 
termined from the overscan region and converted to elec- 
trons by the gain appropriate to each amplifier. A high 
signal-to-noise averaged dark current frame was then sub- 
tracted. Normalized flat fields were derived from exposures 
taken with a halogen lamp immediately after the science 
exposures. For each dither position the average of the other 
two positions was subtracted, to remove the fringing pattern 
(which is approximately fixed in space and does not vary 
significantly with time), while preserving the object signal. 
Residual sky subtraction was performed by fitting a low- 
order polynomial to each detector column (parallel to the 
slit). Individual nod positions were registered and combined 
using a cosmic ray rejection algorithm, and data on the same 
mask from different nights were added using inverse- variance 
weighting. 

One-dimensional spectra were extracted using the IRAF 
apall package: extraction widths of 6 pixels (1.3") were used, 
and the curvature of individual object spectra was "traced" 
across the array. Wavelength calibration was obtained from 
Ne+Ar reference arc lamps, and a fourth-order polynomial 
fit to the centroids of 40 arc lines created a wavelength so- 
lution with rms residuals of 0.3 A. 



3.2 DEIMOS multi-object spectroscopy 

We also obtained slit-mask spectra of 38 EROs using the 
DEIMOS spectrograph (Faber et al. 2003; Phillips et al. 
2002)) at the Nasmyth focus of the 10-m Keck n telescope. 
We used one slitmask each on the CDFS field and the NTT 
field, targeting EROs within the 16.5 x 5 arcmin DEIMOS 
field. Details are given in Table Q Each slit was l"wide, 
with a minimum slit length of 6". DEIMOS has 8 MIT/LL 
2k x 4fc CCDs with 15 /im pixels and an angular scale of 
0.1185" pix -1 . 



Observations were obtained using the Gold 
1200 line mm" 1 grating in first order tilted to a central 
wavelength of 8400 A with a dispersion of 0.32 A pixel -1 . 
Each spectrum spanned w 2600 A and for all targets covered 
AA b s 7000 - 9000 A. A small 8 A region in the centre of the 
wavelength range falls in the gap between two CCDs. We 
used the OG550 order-blocking filter to remove all light at 
wavelengths short-ward of 5500 A, to avoid contamination 
by second-order light. 

The spectral resolution is AAp^ HM ~ 1.4 A 
(Adfwhm ~ 55 km s -1 ), as measured from the sky lines. 
The seeing was typically 0"6-0"9 FWHM, i.e. smaller than 
the slit width, implying that the resolution quoted above is 
a slight overestimate. 

The seeing was typically in the range 0'.'6-0"9 FWHM, 
smaller than the slit width of l'.'O. As the seeing disk was 
smaller than the slit width the true resolution is somewhat 
better for a source which does not fill the slit. 

Observations were divided into individual exposures of 
duration 2400 s. The telescope was dithered 1.5" along the 
slit between integrations. Flux calibration and telluric cor- 
rection were carried out using observations of the spec- 
trophotometric standard HZ 44 (Massey et al. 1988; Massey 
& Gronwall 1990). The flux calibration was checked using 
the spectra of the five alignment stars of known broad-band 
photometry (/ ~ 17 — 19 mag), which were used to position 
the masks through 2" x 2" alignment boxes. Wavelength 
calibration was obtained from Ne+Ar+Hg+Kr reference arc 
lamps. 

Most of the spectra were reduced using vl.1.4 of the 
DEEP2 data pipeline 3 . The pipeline rectifies the slitlets, 
which are curved on the detector, by tracing the slit edges, 
then flat field corrects and determines a 2-D wavelength so- 
lution for each slit. For each of the science frames, a 6-spline 
sky model is fitted and subtracted. The frames are finally 
combined resulting in a mean, sky-subtracted, cosmic ray 
rejected 2D spectrum for each slitlet. 

We found the pipeline failed for around 25% of our 
data, largely because our observational strategy deviated 
from that employed for the more routine DEEP2 survey. 
Our slits were long and untilted; this caused the 6-spline so- 
lutions to diverge. In those cases where the sky subtraction 
from the automated pipeline was inadequate, we re-reduced 
the slitlets using standard IRAF long-slit procedures (as de- 
tailed in Section 13.11 . 



3 http : //astron . berkeley . edu/~cooper/deep/spec2d/ 
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Figure 2. Distribution of EROs for the three LCIRS fields. Can- 

# 1 2te6 t R r J «S l ; t S#NR9?&' l 0©9| : fflflS arc marked with boxes. Where 
redshifts were identified the boxes are marked accordingly (other- 
wise indicated '?'). The areas delineated refer to the spectroscopic 
footprints of the slit-mask. 



4 REDSHIFT DISTRIBUTION AND 
COMPLETENESS ISSUES 

In this section we discuss the determination of spectroscopic 
redshifts and present the redshift distribution for our J — 
H > 3, H < 20.5 sample taking into account various sources 
of incompleteness. 

A complete inventory of all our spectra is given in Fig- 
ures 15, 16, and 17. 



4.1 Redshifts and Spectral Classification 

1-D spectra were extracted over a width of 1.3". Sky resid- 
uals due to imperfect sky subtraction were masked out and 
interpolated over. In both the LRIS and DEIMOS spectra 
typically around ~ 10% of the array was masked out. Al- 
though the DEIMOS resolution is much higher, we used a 
harsher sky threshold in the masking as the background sub- 
traction by the pipeline was not always quite adequate (Sec- 
tion 

To facilitate redshift measurement, we compared with 
early and late type spectral templates from the Gemini Deep 
Deep Survey (GDDS; Abraham et al. 2004) and a lumi- 
nous red galaxy template from the Sloan Digital Sky Sur- 
vey (SDSS; Eisenstein et al. 2003), matching spectral fea- 
tures and the general continuum shape by visual inspection. 
Typically, identifications were based on the 4000A break, 
Call H&K (3968/3933A) absorption, G-band (4300A) ab- 
sorption, or the [O n] AA 3726.2,3728.9 A emission line dou- 
blet which was resolved in the DEIMOS data. 

The redshifts for each of the 3 fields are summarised 
in Tables |H] 31 an d Q Quality flags have been assigned as 
follows: 0=fail, l=uncertain (based on weak feature(s)), 2= 
probable (based on one reliable feature) , 3=secure (based on 
more than one certain feature) and -l=possible high redshift 
(z > 1.4) on the basis of a strong continuum detection but 
no features (discussed below). The redshift distribution for 
the whole sample is shown in Figure |H] Two (5%) of our 
sources are Galactic stars. 

Of 75 objects observed in the three fields, we were able 
to identify redshifts for 44 sources (including two stars) . Of 
these, seven identifications (~ 15%) have the lowest quality 
flag (l=uncertain). Where we could not identify redshifts, it 
was usually because the continuum was too faint to detect 
absorption lines. In five cases we do detect significant contin- 
uum but fail to identify any spectral features (quality=-l). 
Most likely these are higher redshift (z > 1.4) sources where 
the most prominent spectral features are redshifted beyond 
our spectroscopic range. In a few cases (~5%), there were 
data acquisition or reduction problems - such as a source 
falling on top of bad pixel rows on the detector. In SSA22 
two objects fell right at the edge of the slit mask and were 
severely vignetted. In the NTT field, six sources fell on an 
area of the chip which was misaligned due to one of the align- 
ment stars being significantly off centre (the northern most 
part of the mask, Figure |2J and the slitlets therefore missed 
the galaxies. For our statistical analyses we will treat these 
eight (two in SSA22 and six in NTT) sources as unobserved. 
Hence, our sample of observed objects is actually 67 in total, 
for which we obtain 44 redshifts, thereby giving ~ 66% com- 
pleteness. This is comparable to or better than the success 
rates in other spectroscopic surveys of faint galaxies with 
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4.2 Effects of Incompleteness 
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Figure 3. Combined redshift distributions across the three LCIR 
fields for various spectroscopic classes. From bottom to top: all 
galaxies, late type, early types (subdivided in E, E+A (diagonal 
fill) / E+A+e (cross-hatched) / E+e(solid fill)). 



R - K > 5 or I — K > 4 (e.g. equivalent to the K20 sur- 
vey work which has reached a spectroscopic completeness of 
62% to K s < 20, R~K S > 5 (Cimatti et al. 2002, 2003); see 
also surveys described in Yan, Thompson, and Soifer 2004; 
Abraham et al. 2004). We also classified our galaxy spectra 
according to the prominence of diagnostic spectral features, 
as follows: 

i) early- type (E): strong absorption features especially Ca 
II H&K, no detectable [Oil] emission and weak or absent 
Balmer absorption (i.e. H5 with rest-frame equivalent width 
<4A). 

ii) E+A: as early-type (i.) but with enhanced Balmer ab- 
sorption (H<5 with rest-frame equivalent width > 4A (see e.g. 
Van Dokkum & Ellis 2003 ), 

iii) E+e: as early-type but with weaker Ca II H&K and 
moderate to strong [O II] emission, 

iv) late-type: showing prominent [Oil] emission (> 10ct, 
/ ~ 2.0 x 10" 18 erg s" 1 cm 2 ) and no significant Ca II H&K. 

v) E+A+c ('mixed'): resemble the E+A class of galaxy, 
but also exhibit [Oil] emission, i.e. 'mixed' stellar popula- 
tions. 

vi) other: stars and AGN. There is only one AGN in our 
sample, which has been classified by the presence of high 
ionisation lines ([NeV],[NeIII]) 

The split in spectral classifications for those sources 
with redshifts is shown in Table |2] 



A key question is whether the 34% incompleteness correlates 
with I — H colour. For example it might be expected that 
sources with prominent features and large 4000 A breaks 
would yield redshifts more readily. Figure 2] shows how in- 
completeness correlates with colour and magnitude and also 
addresses the question of the extent to which our spectro- 
scopic targets are representative of the overall distribution 
of EROs. The histograms in I — H and H-magnitude show 
numbers of EROs in the detector footprints in all fields, sub- 
divided into those for which we have identified redshifts, not 
identified and not observed. It is clear that the subsample 
of identified redshifts is a representative sample of EROs at 
least to I — H < 3.5 and H < 20, with no bias in colour 
or magnitude space. There is some apparent incompleteness 
towards the redder and fainter end of the distribution, which 
may possibly be populated by z > 1.4 galaxies which are too 
red and faint to be detected optically. 

Of greater concern is the likelihood that some of the 
34% of unidentified sources lie at higher redshifts because 
the principal spectroscopic features are shifted beyond our 
wavelength range, viz. z >1.4. As discussed, a fraction of 
our sample (~ 10%) have quality flags of -1 indicating a 
strong continuum in which features might be expected. A 
key question is whether most or all of these objects lie be- 
yond z ~1.4. 

We tested for this hypothesis in the CDFS field using 
the photometric redshift technique described in Daddi et 
al. (2004). These authors have proposed that galaxies at 
2 > 1.4 will lie in a well-defined region on a BzK colour- 
colour diagram. To facilitate the comparison we took F450W 
(B) and F850LP (z) photometry from the GOODS-South 
HST/ACS (Dickinson et al. 2003; Giavalisco et al. 2004) 
and Ks band data from the ESO ISAAC survey (Vandame 
et al., in prep). The Ks band image covers a 13 of our EROs, 
and of the 6 in this subsample for which we failed to derive 
redshifts, four lie in the BzK region for expected z > 1.4 
(Figure [KJ. However, the region is contaminated by three 
lower z objects, a somewhat higher fraction than that quoted 
by Daddi et al. (2004) . So while this seems to indicate that 
some of our unidentified sources most likely do lie at z > 
1.4, the technique is not definitive for a small sample. The 
result is strengthened however, by the distribution of our 
galaxies in the V — I v's I — H two-colour plane (Figure |BJ, 
which is also generally effective in isolating z > 1.4 galaxies 
(McCarthy 2004). 



4.3 Comparison with Earlier Work 

There have been several redshift surveys of faint objects in 
our target fields, and consequently there is some small over- 
lap (~20%) with our ERO sample. As a check on our identifi- 
cations, we matched our catalogue by coordinates (to within 
1 arcsec), with the GDDS sample for the SSA22 field (Abra- 
ham et al. 2004) and with the FORS2 spectroscopic sample 
of Vanzella et al. (2004) and the K20 survey 5 (Cimatti et al. 
2002) for the CDFS field. For the FORS2 sample (Table EJ, 
we find good agreement with the five sources in common. 



www. arcetri.astro.it/~k20 
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Field 


completeness 


late type 


early type 


E+A 


E+e 


E+A+c 


other 




CDFS 


56% 





45% 


0% 


22% 


22% 


1 star (~ 11%) 




NTT 


87% 


5% 


20% 


20% 


30% 


15% 


1 star (~ 5%), 1 AGN (r 


- 5%) 


SSA22 


54% 


40% 


27% 


13% 


13% 


7% 







Whole Sample 


66% 


16% 


27% 


13.5% 


23% 


13.5% 


2 stars (5%), 1 AGN 


2%) 



Table 2. Classification of sources by field. 





Figure 4. Histograms in (a) I — H and (b) H magnitude showing incompleteness trends in our sample. There appears to be some 
incompleteness towards the red (I — H > 3.5)and faintest (H > 20) end of the distribution, which may be accounted for by possible 
z > 1.4 sources which are too red and/or faint to be detected optically. 



ID 


RA (J2000) 


Dec (J2000) 


I mag 


H mag 


Instrument 


redshift 


quality 


nature 


comments 


939 


3 32 14.98 


-27 42 25.00 


22.45 


18.83 


Iris 











1276 


3 32 19.14 


-27 40 40.40 


22.00 


18.73 


Iris 


1.1290 


3 


E+e 




1430 


3 32 21.15 


-27 41 6.80 


22.40 


19.35 


Iris 


1.1097 


3 


E+A+c 




1447 


3 32 21.24 


-27 40 9.70 


22.55 


19.51 


Iris 


1.0390 


1 


early 




1665 


3 32 24.16 


-27 42 11.00 


21.90 


18.84 


Iris 


0.0000 


3 


M star 




2191 


3 32 30.51 


-27 40 30.40 


23.46 


20.21 


Iris 











2459 


3 32 33.87 


-27 42 4.10 


23.48 


20.42 


Iris 


1.2270 


2 


early 




2588 


3 32 35.62 


-27 43 10.20 


22.74 


19.64 


dcimos 


1.2210 


3 


E+A+e 




2668 


3 32 36.27 


-27 42 49.50 


23.93 


20.18 


dcimos 




-1 




possible z > 1. 


2681 


3 32 37.18 


-27 46 8.20 


21.67 


18.57 


dcimos 


1.0960 


3 


early 




2800 


3 32 39.61 


-27 47 9.00 


23.35 


20.14 


dcimos 


1.3180 


2 


early 




3695 


3 32 41.64 


-27 41 52.10 


23.09 


19.92 


dcimos 











3807 


3 32 43.92 


-27 42 32.40 


23.67 


19.67 


dcimos 




-1 




possible z > 1. 


4071 


3 32 42.81 


-27 43 28.90 


23.80 


20.10 


dcimos 




-1 




possible z > 1. 


1792 


3 32 25.77 


-27 43 47.60 


23.48 


20.28 


dcimos 











2158 


3 32 30.33 


-27 45 23.60 


22.54 


19.44 


dcimos 


1.2220 


3 


E+e 





I1843 4 


3 32 26.28 


-27 45 36.00 


22.78 


19.93 


deimos 


1.2250 


3 


E+A 


12897 


3 32 38.12 


-27 44 32.90 


22.48 


19.65 


dcimos 


1.2210 


2 


early 


12331 


3 32 32.97 


-27 41 16.90 


21.66 


18.68 


dcimos 


1.0420 


3 


E+A+c 


13063 


3 32 39.51 


-27 41 17.30 


22.71 


19.83 


dcimos 


1.0394 


3 


early 


12876 


3 32 38.43 


-27 40 19.40 


22.72 


19.90 


dcimos 


1.0370 


3 


early 


12163 


3 32 30.00 


-27 47 26.40 


23.62 


20.19 


dcimos 


0.9900 


1 


early 



Table 3. CDFS Sample: Quality flags: 0=fail, l=uncertain, 2=probable, 3=secure, -l=z> 1.4,6=star. Objects prefixed with an T in 
the above table, were selected with R — K > 5. 
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ID RA (J2000) Dec (J2000) I mag H mag I — H redshift quality nature 



275 


12 


04 


22.73 


-07 


31 


15.6 


22 


.90 


19.96 


2.94 






-1 


poss hi-z 


716 


12 


04 


34.60 


-07 


30 


46.3 


23 


.05 


19.99 


3.06 


1 


.1865 


3 


E+A 


1193 


12 


04 


35.45 


-07 


30 


18.9 


22 


.78 


19.69 


3.09 


1 


.2680 


3 


E+e 


1218 


12 


04 


32.04 


-07 


30 


10.9 


22 


.56 


19.55 


3.01 






-1 


poss hi-z 


1388 


12 


04 


19.58 


-07 


29 


54.2 


23 


.59 


20.07 


3.52 


1 


.2700 


3 


E+e 


1494 


12 


04 


26.70 


-07 


29 


46.3 


23 


.13 


19.84 


3.29 


1 


.2680 


3 


E+e 


1496 


12 


04 


18.98 


-07 


29 


45.4 


23 


.45 


20.50 


2.95 


1 


.1020 


2 


early 


1962 


12 


04 


24.95 


-07 


29 


07.4 


23 


.41 


20.21 


3.20 






-1 


poss lii-z 


2159 


12 


04 


28.38 


-07 


28 


53.0 


23 


.44 


20.07 


3.37 


1 


.2680 


3 


E+A 


3249 


12 


04 


14.99 


-07 


27 


43.4 


22 


.41 


19.37 


3.04 


1 


.1000 


3 


E+A 


3300 


12 


04 


19.11 


-07 


27 


33.4 


23 


.59 


20.31 


3.28 


1 


.0660 


3 


early 


3456 


12 


04 


18.77 


-07 


27 


27.2 


22 


.71 


19.51 


3.20 


1 


.0700 


3 


E+A+c 


3464 


12 


04 


08.21 


-07 


27 


26.3 


22 


.27 


19.35 


2.92 


1 


.1750 


3 


AGN(?) 


3472 


12 


04 


18.22 


-07 


27 


18.1 


23 


.39 


19.96 


3.43 


1 


.0740 


3 


E+A 


3720 


12 


04 


16.12 


-07 


27 


01.2 


23 


.49 


20.15 


3.34 


1 


.0700 


3 


early 


4206 


12 


04 


17.36 


-07 


26 


26.2 


22 


.69 


19.36 


3.33 


1 


.1450 


3 


E+e 


5928 


12 


04 


02.05 


-07 


21 


18.9 


23 


.06 


20.15 


2.91 


1 


.3185 


3 


E+A+c 


5962 


12 


04 


07.80 


-07 


24 


15.9 


23 


.04 


20.12 


2.92 


1 


.1000 


3 


E+e 


6052 


12 


03 


57.36 


-07 


24 


15.9 


23 


.06 


20.14 


2.92 





.9600 


1 


early 


6056 


12 


04 


01.39 


-07 


24 


09.6 


23 


.00 


20.01 


2.99 


1 


.2230 


3 


E+A+c 


6322 


12 


03 


56.88 


-07 


23 


58.9 


22 


.24 


19.20 


3.04 


1 


.1030 


3 


E+e 


6696 


12 


03 


55.82 


-07 


23 


30.1 


22 


.64 


19.27 


3.37 


1 


.2180 


1 


late 


7463 


12 


03 


45.06 


-07 


22 


32.0 


23 


.71 


18.99 


4.72 





.0000 


6 


star 



Table 4. NTT Sample 



ID 


RA (J2000) 


Dec (J2000) 


I mag 


H mag 


mask# 


redshift 


quality 


nature 


403 


22 


17 


49.01 


22 14.30 


23 


.46 


20.44 


2 









527 


22 


17 


41.50 


21 43.60 


23 


.58 


20.29 


2 









531 


22 


17 


44.84 


21 42.00 


23 


.58 


20.40 


2 









713 


22 


17 


40.86 


20 54.50 


23 


.51 


20.42 


2 


1.2890 


2 


late 


749 


22 


17 


44.03 


20 45.90 


23 


.52 


20.47 


2 









781 


22 


17 


47.46 


20 26.90 


22 


.04 


18.86 


2 


0.9146 


3 


E+e 


958 


22 


17 


46.73 


19 53.90 


23 


.50 


20.49 


2 


1.2252 


3 


E+A+c 


1051 


22 


17 


48.87 


19 26.90 


21 


.80 


18.64 


2 


1.2280 


2 


early 


1267 


22 


17 


49.36 


18 39.60 


23 


.26 


20.19 


2 


1.1076 


3 


E+A 


1535 


22 


17 


35.46 


17 34.60 


23 


.45 


19.81 


1 









1623 


22 


17 


48.27 


17 12.30 


22 


.92 


19.86 


1 









1655 


22 


17 


48.59 


17 5.80 


23 


.39 


19.85 


2 


1.1990 


1 


early 


1792 


22 


17 


45.73 


16 36.60 


23 


.31 


19.55 


1 









1803 


22 


17 


34.62 


16 36.30 


23 


.37 


20.30 


1 









1920 


22 


17 


48.37 


16 9.00 


23 


.69 


20.10 


2 


1.1420 


1 


late 


2220 


22 


17 


41.04 


15 1.30 


23 


.89 


20.30 


2 


0.7550 


1 


late 


2244 


22 


17 


43.02 


14 51.00 


23 


.46 


19.77 


1 









2260 


22 


17 


32.86 


14 48.10 


22 


.99 


19.50 


1 


0.7980 


2 


early 


2335 


22 


17 


42.08 


14 30.80 


23 


.37 


19.91 


1 









2467 


22 


17 


35.79 


13 52.40 


22 


.12 


18.99 


1 


1.0320 


3 


E+e 


2593 


22 


17 


47.54 


13 27.50 


23 


.01 


19.34 


1 









2638 


22 


17 


41.71 


13 18.00 


23 


.32 


20.21 


1 


0.7649 


3 


E+A 


2767 


22 


17 


39.50 


10 46.50 


22 


.78 


19.56 


1 









2779 


22 


17 


47.18 


10 41.30 


22 


.89 


19.87 


1 


1.0220 


2 


late 


2993 


22 


17 


44.86 


12 6.70 


22 


.91 


19.88 


1 


1.3695 


2 


late 


3116 


22 


17 


43.25 


13 1.10 


23 


.23 


19.82 


1 









3253 


22 


17 


32.34 


12 52.40 


23 


.25 


20.15 


1 


1.1530 


2 


late 


3378 


22 


17 


33.67 


13 5.30 


23 


.38 


19.54 


1 


0.7620 


1 


early? 



Table 5. SSA22 Sample 
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galaxies at z>1.4 / Old galaxies at z>1.4 




(B-z) 



Figure 5. Colour-colour (z — K)ab vs (B — z)ab diagram for 
galaxies in the CDF-S. According to Daddi et al. (2004), the 
region defined by (z — K)ab — (B — z)ab > —0.2 (blue shaded 
area) should isolate z > 1.4 star-forming galaxies. Conversely, 
the intersection of the regions (z — K)ab — (B — z)ab < —0.2 
and (z — K)ab > 2.5 should isolate old galaxies with z > 1.4 
(yellow shaded area). 4 out of 6 galaxies without rcdshifts (marked 
with an asterisk) lie in these regions (although the contamination 
from lower redshift sources is somewhat higher than the fraction 
quoted by Daddi et al. (2004). Sources flagged as potential z > 1.4 
sources are overlaid with boxes - 2 out of 3 fulfill the necessary 
colour criteria. 




3.0 3.5 4.0 

I-H (LCIRS,Vega) 



4.5 



Figure 6. V — I v's I — H colour-colour diagram. Only galax- 
ies with spectroscopically identified redshifts with class > 1 are 
shown. Open diamonds are z < 1 galaxies, filled diamonds are 
1 < z < 1.4, and large filled circles are unidentified but believed 
to be at z > 1.4 on the basis of featureless continuum. The remain- 
ing unknown redshifts are marked '?' and stars are marked '*'. 
Galaxy tracks are marked to give an indication of where z > 1.4 
galaxies might be expected to lie. The tracks are dotted at 2 < 1.3 
and solid at z > 1.3. 



Four have high confidence quality flags in both samples and 
agree to within Az = 0.001. The redshift for the remaining 
source is unidentified in both surveys. Of the four sources 
overlapping with the K20 survey (Tabled, three of the red- 
shifts are in agreement, but one which we have flagged as low 
confidence at z—0.99 is discrepant with the higher z=1.553 
K20 estimate. 

Matches with the GDDS sample are shown in Table |S] 
Of the three overlapping sources for which we have confi- 
dent identifications, one agrees and two disagree. The GDDS 
confidence rating for one of these failures is < 50%. In three 
cases where our redshift identification is insecure (class 1) 
and in two cases where we failed, GDDS find higher redshifts 
(2 > 1.395) with a high confidence rating, supporting our 
hypothesis that the redshifts which we failed to get poten- 
tially lie at higher- z . 



5 ANALYSIS 

We now turn to addressing the key objectives for the survey 
set out in ^ Summarizing the situation, using the photo- 
metric criteria H < 20.5 and / — H > 3.0, we have drawn 
spectroscopic targets from 252 EROs in three LCIR sur- 
vey fields, of which 159 fall within the field of view of the 
slit masks employed. We successfully observed a sub-sample 
of 67 objects drawn randomly from the photometric sample. 
Spectroscopic identifications were secured for 44 sources (in- 
cluding two stars), a completeness fraction of 66%. We have 
demonstrated that this large spectroscopic sample is repre- 
sentative of the ERO population in the LCIR survey. The 
spectra are generally of good quality with higher dispersion 
than in previous studies (e.g. K20, GDDS). 



5.1 Redshift Distribution 

Our redshift distribution (N(z), Figure [3J ranges between 
0.755 and 1.4 (with some sources potentially at 2 > 1.4), 
and therefore demonstrably supports the hypothesis that an 
/ — H > 3 colour criterion is effective in selecting galaxies 
at z > 1 (McCarthy et al. 2001). 

Figure |7] compares the observed distribution and that 
predicted for the larger LCIRS sample using photometric 
redshifts based on UBVRIH photometry in the indepen- 
dent LCIRS field HDF-S (Firth et al. 2002). From the lat- 
ter distribution, Firth et al. concluded that popular semi- 
analytic models provided a reasonable fit to the overall N(z) 
but under-predicted the abundance. 

The redshift distribution in our spectroscopic survey 
varies from field to field with clear structures in each, sug- 
gesting significant cosmic variance and clustering among 
the population. Furthermore, different mixtures of spectral 
types are seen in the various fields (Table |2Jl. The simplest 
conclusion consistent with this disparity is that the domi- 
nant component of the ERO population is made of luminous 
spheroidals which are strongly clustered. We will return to 
analyzing the variance seen across the three fields when we 
measure the overall spatial clustering in our sample in § 15.51 

Although our spectroscopic redshift distribution is lack- 
ing the high— z tail seen in the photometric distribution - 
due to our incompleteness in the 'spectroscopic desert' - or 
2 > 1.4 regime, we find an overall abundance in agreement 
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ID 


redshift 


I mag 


FORS2 ID 


FORS2 redshift 


FORS2 I mag 


1276 


1.129 


22.00 


J033219. 15-274040.2 


1.128 


22.25 


2681 


1.096 


21.67 


J033237. 19-274608.1 


1.096 


21.98 


2800 


1.318 


23.35 


J033239.64-274709.1 


1.317 


23.69 


2158 


1.222 


22.54 


J033230.34-274523.6 


1.223 


22.89 


1792 




23.48 


J033225. 76-274347.0 




24.13 



Table 6. Redshift comparison with the FORS2 survey in the CDFS field. 



ID 


redshift 


quality flag 


K20 id 


K20 redshift 


2681 


1.09600 


3 


CDFS.00633 


1.09600 


2158 


1.22200 


3 


CDFS.00507 


1.22300 


11843 


1.22500 


3 


CDFS.00547 


1.22200 


12163 


0.990000 


1 


CDFS.00139 


1.55300 



Table 7. Redshift comparison with the K20 survey in the CDFS field. 



3.0 - 
2.5 - 




redshift 



Figure 7. Observed redshift distribution for 3 LCIRS fields, cor- 
rected for incompleteness (histogram - with shaded region pos- 
sible z > 1.4 galaxies), compared with that estimated photo- 
metrically by Firth et al. (2002; dashed curve) in a fourth field. 
Overplotted is a predicted N(z) if EROs produced all ellipticals 
today via passive luminosity evolution. The local elliptical pop- 
ulation has been reverse evolved (assuming all are 5Gyr old), 
using the luminosity function derived by Nakamura et al (2003) 
for E/S0 galaxies in the SDSS, which has a = —0.83 (slope), 
0*=O.OO47*(O.7) 3 per Mpc 3 , M*= -20.75+5.*alogl0(0.7)) The 
solid curve is the backwards-PLE model including our I — H > 3.0 
colour cut, and dotted curve has no colour cut on PLE). 

with that found by Firth et al. in the HDF-S, i.e. signifi- 
cantly more than is predicted by semi-analytic models. 

5.2 I — H versus R — K Colour Selection 

Our / — H > 3 selection criterion is different from the 
frequently-used R-K > 5 (Cimatti et al. 2002) or R-K > 6 
(Hu & Ridgway 1994) cut and thus it is interesting to con- 
sider whether the greater prominence of early-type spectra 
(and the larger associated cosmic variance) is due, in part, to 



this difference. Our fraction of actively star-forming galaxies 
(~ 16%) is much less than that found (50%) in the R — K 
selected K20 sample (Cimatti et al. 2002) or the bright K- 
sample of Yan et al. (2004) . 

We also might reasonably expect the shorter wave- 
length baseline of I — H (compared to R — K) to provide 
an increased sensitivity to strong age-dependent continuum 
breaks as opposed to the broader wavelength signatures aris- 
ing from reddening. Unfortunately, given the different na- 
ture and classification methods for the various spectroscopic 
datasets, it is unclear whether the different fractions of early- 
type galaxies arises entirely from the colour selection. Part of 
the difference may arise from the increased signal to noise 
in our data. For example, although a further 40% of our 
galaxies reveal modest [O n] emission, our spectra enable us 
to clearly demonstrate an underlying old stellar population, 
whereas such galaxy types may previously have been classi- 
fied as dust-reddened starbursts. Figure 1 also demonstrates 
that E+A galaxies, such as appear in our sample (for exam- 
ple, with 5% burst mass in the last 100 Myr) would not be 
expected in a selection based on R — K > 5 colour cut. 



5.3 Spectroscopic Properties of the ERO 
Population 

We now turn to the important question of the nature of 
the ERO population as revealed by the Keck spectra. Thus 
far we have categorized each spectrum according to its ap- 
pearance, deferring any astrophysical interpretation of the 
various classifications. Our goal is to address the relevance 
of the E, E+A, E+e, E+A+e and late types introduced in 
§4.1 as well as the homogeneity of the overall ERO popu- 
lation. We stress that our analysis in this section will be 
illustrative given that the spectroscopic data alone cannot 
separate uniquely many of the key variables e.g. the effects 
of age, metallicity and complex star formation histories (c.f. 
McCarthy et al. 2004). 

We begin by stacking the spectra according to our ba- 
sic classification scheme. Figure |3] shows the result of this 
coaddition for the higher quality spectra selected with qual- 
ity flags 2 or 3. In practice, given we used both DEIMOS 
and LRIS spectrographs, we restricted the coaddition of the 
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ID 


redshift 


quality flag 


GDDS ID 


GDDS redshift 


GDDS I mag 


GDDS confidence 




3253 


1.153 


2 


23.25 


SA22-0062 


1.154 


23.17 


[O ii] 


2260 


0.798 


2 


22.99 


SA22-0448 


1.202 


22.80 


<50% 


2467 


1.032 


3 


22.12 


SA22-2548 


1.022 


21.96 


certain 


1655 


1.199 


1 


23.39 


SA22-0674 


1.493 


23.19 


>75% 


1920 


1.142 


1 


23.69 


SA22-1983 


1.488 


23.68 


95% 


3378 


0.762 


1 


23.38 


SA22-0107 


1.448 


23.20 


single emission line 


2335 







23.37 


SA22-0398 


1.395 


23.16 


95% 


2593 







23.01 


SA22-0189 


1.490 


22.80 


95% 


1535 







23.45 


SA22-0948 


1.396 


23.24 


none/best guess 


1792 







23.31 


SA22-0721 


1.483 


23.14 


none/best guess 


1803 







23.37 


SA22-0717 


2.060 


23.10 


none/best guess 



Table 8. Redshift comparison with the GDDS survey in the SSA22 field. 




3800 3900 4000 4100 4200 




no t Li 

3800 3900 4000 4100 4200 

X m , I Ang 



Figure 8. Noise spectrum for the E (top) and E+A (bottom) 
composite spectra, composed of the individual spectra normalised 
to the same continuum flux and overlaid. The standard deviation 
in this spectrum is then the error on the composite spectrum. 



E, E+A and E+e spectra, for which absorption line mea- 
sures are particularly helpful, to those taken with the supe- 
rior resolution DEIMOS, rebinning to a final dispersion of 
S.SApixeP 1 . In the case of the late-type spectra, as most 
examples occur in the LRIS sample, we restricted the coad- 
dition to those taken with LRIS rebinned to 13.6Apixel _1 . 
We find the error on each composite spectrum to be the de- 
viation in the spectra when each is normalised to the same 
relative flux in the wavelength range 3800 — 420CL4. The noise 
spectra for the E and E+A galaxies are shown in Figure [S] 
we use this to deduce the error on diagnostics such as the 
D4000 break and line equivalent widths. 

We now briefly discuss the coadded spectra in turn: 



5.3.1 Pure early type (E) 

About 30% of our sample are in the E category and our 
luminosity-weighted mean spectrum has a 4000A break 
D 4000 ~ 1.7 6 . Comparing this with the hypothesis of a sim- 
ple stellar population following a single burst in the context 
of Bruzual & Chariot (2003) models we find an average age 
of ~ 2 Gyr at the mean redshift of ~z—1.2, assuming solar 
metallicity (Figure HTk ) . This implies a formation epoch of 
2 > 2 for this subset of the ERO population. This result 
agrees with McCarthy et al. (2004) who find a < z f >=2.4 
for 20 red galaxies with z > 1.3 in the GDDS, and deduce 
an early and rapid formation for a substantial fraction of 
these. Furthermore, Cimatti et al. (2002,2003,2004) deduce 
formation redshifts between 2-3.4 for similar objects found 
in the K20 survey. 

5.3.2 Intermediate age (E+A) 

A significant fraction of our sample have strong Balmer ab- 
sorption lines in addition to prominent 4000A breaks and are 
indicative of secondary star-formation. These galaxies have 
traditionally been referred to as E+A galaxies (also known 
as 'K+A'; Dressier et al. 1999 ; see also Blake et al. 2004 
and refs. therein). Although ~ 30% of our absorption line 
sample fall in this category, the fraction shows considerable 
field-to-field dispersion (0-50% across the 3 LCIRS fields). 
Nonetheless, the rate of occurrence is much higher than in 
nearby field galaxies. 

Zabludoff et al. (1996) find 0.6% of all galaxies at z ~ 
in the Las Campanas Redshift Survey to be E+A (as se- 
lected by EW(H<5)> 4.5A). Given that EROs make up about 
~ 10% of the overall galaxy population at our limiting mag- 
nitude, we find more than twice as many E+As at z ~ 1. 
With a stricter criterion, (EW(HJ) > 5 A) Goto (2005) finds 
a local E+A incidence of ~ 0.1% in the Sloan Digital Sky 
Survey (SDSS; York et al. 2000). Bearing in mind that our 
selection function for EROs is picking out only E+A galax- 
ies which are spheroidal in origin (locally, E+As are found 
to also have late-type spiral progenitors, see e.g. Blake et al. 
2004, Zabludoff et al. 1996), this is reasonably strong evi- 
dence for evolution in the number density of such systems 
from z~l-2toz~0. 

6 We define the index D4000 = /(3750-3950A)//(4050-4250A) 
((Bruzual 1983) 
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Figure 9. Stacked spectra for categories E, E+A, E+e and late type as denned in §4.1. The E, E+A and 
the superior resolution DEIMOS data whereas the late-type spectra were restricted to the LRIS sample; Only specte 
2 or 3 were included. 
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Figure 10. (a) Pure early type composite spectrum compared to a Bruzual and Chariot (2003) model spectrum (dashed line) of a 
3.5Gyr population and (b) the E+A spectrum is well fit by an underlying 2Gyr population, with a lOOMyr old secondary burst of 2% 
by mass. Both assume a Salpeter IMF and solar metallicity. The shaded regions show the errors on the composite spectra, which are the 
standard deviation of the flux at each wavelength in the individual spectra (scaled to the same relative luminosity). 



The composite spectrum (Figure |S| has an average 
equivalent width Whs = 6A. Figure lllb shows the evolu- 
tion with age of the H<5 equivalent width for sub-solar, solar 
and super-solar metallicities. Combining with the D4000 in- 
dex (S I5.3.1I Figure lT11 we find the spectrum can be readily 
interpreted in terms of an underlying population as old as 
that of the mean E type spectrum discussed above, with an 
additional secondary component (5-10%) which formed be- 
tween 30 Myr and 1 Gyr prior to the epoch of observation. 
This is consistent with examples of star formation activity 
found in 2 ~ 1 spheroidals by Treu (2004 and references 
therein). 



In summary, the simplest hypothesis is that the E+A 
galaxies have a similar origin to the E-type galaxies reflect- 
ing that component of the population which recently under- 
went a secondary burst, possibly associated with a merger. 
It is possible that the K-bright star forming objects at 2 ~ 2 
seen in GOODS and K20 could be the same objects we are 
seeing as E+As at 2 ~ 1, in their post-burst phase. 



5.3.3 Active early type (E+e) 

As earlier studies in clusters have shown (e.g. Couch & 
Sharpies 1987, Barger et al 1996), for established stellar 
systems undergoing secondary activity, a subset are likely 
to show the mixed signals of an old stellar population and 
emission lines during the active phase. In the case of merger 
driven evolution, unless the accretion is purely stellar such 
systems will show diluted evidence of an old stellar pop- 
ulation and [O II] emission. Depending on the duty-cycle, 
the Balmer absorption lines will be partially or wholly filled 
by gaseous emission. In fact, although Balmer absorption 
lines are undetected amongst each indivual spectrum in this 
class, they are weakly visible in the composite spectrum 
(Figure I^J. 

Given our explanation for the E+A sources, it is 
natural to adopt the above explanation for many of the 
E+e spectra whose line properties are understandably more 
diverse. 

In summary, the dominant component of our ERO sam- 
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Figure 11. (a) Evolution of D4000 according to a single- burst Bruzual & Chariot (2003) model of solar metallicity (solid line), 0. 2 X solar 
(dotted line) and 2.5x solar (dashed line). Comparison with the mean early-type spectrum implies an age of ~ 1.6 Gyr at the mean 
sample redshift 2 ~1.2. (b) Evolution of EW(H<5) for a single-burst Bruzual & Chariot (2003) model of solar metallicity (solid line), 
0.2xsolar (dotted line) and 2.5x solar (dashed line). 



pie consists of old galaxies seen in various stages of inter- 
mittent, but minor, activity. Figure l"i"2l shows an exempla- 
tive star formation cycle which may give rise to the var- 
ious stages of evolution we observe in our sample. The E 
type spectra are old galaxies are formed at high redshift, 
perhaps in an initial discrete burst of star formation. Some 
time later a secondary burst is triggered, giving rise to the 
[Oil] emission seen in the E+e spectra, and subsequently 
Balmer aborption lines become visible once star formation 
has been truncated, in the E+A class. The mixed popula- 
tions seen, which are predominantly E+A with [Oil] emis- 
sion, represent some intermediate stage between E+e and 
E+A. Barger et al. (1996) successfully explain the numbers 
of systems similar to those discussed above in clusters at low 
redshift (z=0.31) using models involving secondary bursts of 
star formation over the 2 Gyr prior to the epoch of obser- 
vation. Treu et al. (2005) have recently demonstrated that 
the star formation history of field spheroidals is heavily de- 
pendant on their mass. If galaxy formation is governed by 
a 'down-sizing' mechanism (Cowie et al. 1996), i.e. growth 
through rapid star-formation trends smoothly to less mas- 
sive galaxies with decreasing redshift, then in a large sample 
we might expect the E, E+A, E+A+c and E+e galaxies to 
be progressively ordered in decreasing luminosity (e.g. in the 
rest-frame /—band). We examined our data, unsuccessfully, 
for evidence of such a trend. It may be that all of these 
galaxies formed the bulk of their mass at high-z and that 
subsequently environment plays a bigger role, with mergers 
or interactions triggering more recent complex star forma- 
tion histories. 

5.3.4 Late T yve 

A small fraction (~ 16%) of our sample are classed as late 
type on the basis of prominent [O 11] emission together with 
weak or non-existent Ca II H&K absorption. In these sys- 
tems, the luminosity weighted stellar population is clearly 
young and star-forming. 

In the co-added LRIS spectrum of 4 late- type galaxies in 



SSA22, we find an average line flux of 1 x 10~ 17 ergs -1 cm~ 2 . 
At the average redshift z=1.2, this equates to a fairly mod- 
est star formation rate (SFR) of 1.1 ± 0.3 MQyr -1 using 
Kennicutt's (1998) conversion: 

SFRiMoyr' 1 ) = (1.4 + 0.4) x KT 41 L ([Oil]) {ergs- 1 ) 
However this may be a lower limit as it is uncorrected 
for extinction which presumably is responsible for the red 
I — H colour. As an illustration, a 'pure' dusty starburst 
(with no significant flux contribution from an underlying 
older population) our I — H > 3 selection criterion requires 
E(B — V) > 1 (Figure 0. In the simplest case where dust 
screens the gas and stars by equal amounts, this would imply 
a correction to the SFR as large as > 50 (~ 100 Moyr -1 ). 
in practice, of course, there are numerous uncertainties in 
deriving star formation rates in this manner. The recent 
Spitzer results of Yan et al. (2004) imply high extinction for 
the late type red galaxies and typical star formation rates 
50 - 17OM yr~l. 



5.4 Abundance of Galaxies with Established 
Stellar Populations 

The foregoing discussion has revealed that the bulk of our 
I — H >3 sample, specifically most of those in the dom- 
inant spectral categories E, E+A and E+e, contain well- 
established stellar populations (defined broadly as those 
where the bulk of the star formation occurred prior to z=2). 
It is thus interesting to make a rough comparison of the lumi- 
nosity density in this population with that in the population 
of present day spheroidals. 

With a surface area of 200 arcmin -2 across our 3 LCIRS 
fields, and assuming a uniform redshift-dependent selection 
function across the interval 1.0< z <1.35 (in which 90% of 
our objects lie) we effectively survey a co-moving volume of 
~ 1.70x 10 4 Mpc 3 . The average space density of EROs is thus 
(f> ~ 1 x 10 -3 Mpc -3 of which we estimate ~ 75% contain an 
established stellar population. At the median redshift z—1.2, 
our H <20.5 magnitude limit corresponds to an absolute 
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Figure 12. Diagram showing a cycle of star formation activity which would explain the observed spectra of the E, E+e, E+A and E+A 
with [Oil]. The E type galaxies are formed in an initial discrete burst of star formation, some time later a secondary burst gives rise to 
the [Oil] emission seen in the E+e spectra and subsequently Balmer aborption lines once star formation is truncated. 



rest-frame I-band luminosity of Mr =-21.85, equivalent in 
local terms (Blanton et al. 2003 ) to ~1.25L*. 

According to the Bruzual & Chariot (2003) models dis- 
cussed in §5.3, an established stellar population will pas- 
sively fade by ~1 magnitude in rest-frame I since z ~1.2 
(e.g. Fontana et al. 2004). Accordingly, in the absence of 
growth by accretion (Treu et al 2004), our survey sensitiv- 
ity corresponds to locating the progenitors of current epoch 
0.6L* galaxies. Our inferred abundance of <f)(L > 0.6L*) = 
0.6(f)* — 3 x 10~ 3 Mpc -3 , or just over a half of all galaxies 
today. 

We thus find agreement with Firth et al. (2002) that the 
ERO population is sufficient to explain ~ 50% the present- 
day spheroidal population solely via passive evolution. This 
is consistant with purely photometric determinations of the 
evolving red rest-frame luminosity function (e.g. Chen et al. 
2004) and both spectroscopic and photometric measures of 
the evolving mass density in red systems (e.g. Glazebrook 
et al. 2004; Bell et al. 2004; Fontana et al. 2004). While the 
most up-to-date semi-analytic hierarchical merging models 
(e.g. Sommerville et al. 2004) can produce the required num- 
ber density of massive systems at z ~ 2, they have great 
difficulty in reproducing the high space density of massive 
galaxies with red colours and evolved stellar populations. 

5.5 Clustering of the ERO Population 

One of the motivations for defining the redshift distribution 
of the ERO population is to better understand their spa- 
tial clustering. McCarthy et al (2001) measured the angular 
clustering of the I — H >3 population in the HDF-S and 
CDF-S field, a total area of 0.39 deg 2 . Within the range ap- 
propriate for verification via our study (19< H <20.5), they 
found an amplitude 8o = 6.7 ± 0.4 arcsec for the two-point 
correlation function (assuming a form w(6) = (6q / 6) 7 " 1 
where the slope of the correlation function 7 is -1.8 as in 
local samples). To invert the angular function and its nor- 
malization 80 into its spatial equivalent for comparison with 
local sources, McCarthy et al adopted a Gaussian N(z) with 
a peak at 2=1.2 and a dispersion <r z =0.3 . They obtained a 
clustering length ro=9.8 Mpc - a value considered to be 
larger than that for present day galaxies but comparable to 
that for the local elliptical population (Norberg et al 2001). 

Firth et al (2002) subsequently analyzed the angular 
clustering more rigorously in the HDS-S field. For a H <20.5 
sample of 170 galaxies with I — H > 3 across a smaller area 



of 0.21 deg 2 a more refined photometric redshift distribution 
was derived. They relate the correlation length ro at a given 
redshift z to its zero redshift equivalent, say r z o, via the 
expression: 

r„( 2 ) = r z0 (l + z)^~ 3)h 

where e parameterizes the evolution of clustering. e=0 
corresponds to constant clustering in proper space and e = 
-1.2 constant clustering in comoving coordinates. Assuming 
their ERO population did not evolve in comoving terms (i.e. 
e=-1.2), Firth et al deduced the population would have a lo- 
cal clustering with r z o=7.5 ±3.7, i.e. 1.5 times larger than 
local galaxies but similar to that observed (r z o ~6-8 Mpc) 
for local ellipticals. Most significant of all, the clustering of 
the ERO sample in both studies is a factor of 2.5-4 times 
stronger than the full H-selected population. The larger er- 
ror bar on the Firth et al determination compared to Mc- 
Carthy et al is only partly due to the smaller field employed; 
it also reflects a more realistic measure of the overall inver- 
sion uncertainty. 

The availability of spectroscopic redshifts allows us to 
strengthen these conclusions in two ways. Firstly, the im- 
proved precision of a spectroscopic redshift over its photo- 
metric equivalent permits us to directly detect large scale 
structures composed of EROs rather than relying entirely on 
projected angular statistics. Secondly, the deduced redshift 
distribution, N(z) (see §4) allows us to significantly improve 
upon the clustering analyses undertaken by McCarthy et al 
and Firth et al. 

To illustrate the former, our spectroscopy, when com- 
bined with those selected with 7? — K >5 in the FORS2 
sample (Vanzella et al (2004) ) define a prominent over- 
density of EROs at 2=1.22 in the CDFS (Figure ED- The 
angular distribution of 20 sources within 8z = ±0.01 (i.e. 
a velocity dispersion of ~ 400 km s _1 ) is shown in Fig- 
ure EH an d appears to trace a large wall- like structure pos- 
sibly consistent with an assembling cluster. The five LCIRS 
galaxies in the structure all exhibit signs of a prominent old 
stellar population, whereas the majority in the FORS2 sam- 
ple were spectroscopically identified by [O 11] emission alone. 
This discrepancy may be explained by our selection function 
for EROs. However, assuming that this is a gravitationally 
bound cluster, we speculate that the relative spatial concen- 
trations of the two galaxy types might also be evidence for 
the morphology-density relation (Dressier et al. 1997 ) oper- 
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Figure 13. Large Scale Structure in the CDFS (including sources 
selected using R—K > 5). The bin size is Az = 0.02 (400km s _1 ). 
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Figure 14. Angular distribution of galaxies around the z = 
1.22 ± 0.01 structure in the CDFS. Each galaxy is labelled with 
its spectroscopic redshift. Crosses represent sources in the FORS2 
sample, diamonds those from the present survey. 

ating at higher redshift. In summary, we are possibly seeing 
the older, evolved galaxies at the centre of the cluster and 
the more active star-forming galaxies tracing the outskirts. 
Assuming further that the cluster is virialized, we find a 
mass of ~ 5.6 x 1O 14 M0 within a radius of 5 Mpc. This is 
consistent with other high redshift clusters , for example, the 
rich cluster MS1054-03 which has a mass of ~ 1.9 x 1O 15 M 
within ~ 2 Mpc (Tran et al. 1999). 

Returning finally to the spatial correlation function, we 
can verify more precisely the discussion in Firth et al (2002) 
by comparing the photometrically-inferred redshift distri- 
bution N(z) for a H < 20.5, I — H >3.0 sample (Figure 
29 in Firth et al's paper) with our spectroscopic equivalent 
(Figure 6). As we have already noted, the distributions are 
remarkably similar so we can expect little change in the 
analysis of Firth et al. 

However, taking the angular clustering signal from the 
larger area studied by McCarthy et al (2001), we infer an 
amplitude for the angular correlation function w(9) at #=1.0 
arcmin of 4=0.173 ±0.010 (c.f. 0.20 ± 0.06 from Firth et 



al's HDF-S study) and using the redshift distribution of Fig- 
ure 6, this inverts to a local clustering scale length for the 
comoving case of r z o= 6.5 ± 0.4. This is somewhat smaller 
than the amplitude claimed by Firth et al but a more precise 
estimate. 

Our measurement is slightly smaller, but potentially 
more accurate, than most ro amplitudes in the literature re- 
lying on photometric redshift estimates. For example, Daddi 
et al. (2002) quote a range 5.5 £ r /(/i _1 Mpc) < 16 (note 
they use Ho = 100/t) for old passively evolving EROs. Brown 
et al. (2005) estimate the spatial clustering of EROs over 
a large area (0.98 2 degrees) in the NOAO Deep Wide Field 
Survey, using photometric redshifts, and extrapolate a value 
of r ~ 7.5/i _1 Mpc for K > 20, where H = 100ft. 



6 CONCLUSIONS 

We have used J/-band imaging from the LCIRS survey to 
select a large and uniform sample of EROs brighter than 
H > 20.5 with a colour selection (I — H) > 3, over three 
fields (SSA22, CDFS & NTT Deep Field). 

Our deep Keck LRIS & DEIMOS spectroscopy has 
targetted 50% of these EROs over a sampling area of 
200 arcmin 2 . Of the 67 EROs appearing on our slitmasks, 
we have determined reliable spectroscopic redshifts for 44 
- a completeness fraction of 66%. This is the most exten- 
sive spectroscopic study of this population to date, and the 
first to focus on an i/-band selection. Most of our spectra 
have a continuum signal/noise of 2-3 per A and the superior 
resolution of our DEIMOS data enables us to undertake di- 
agnostic spectroscopic studies of various subsets of the ERO 
population. 

We find that most of the ERO population contain a 
clear signature of an old stellar population. In contrast to 
other studies, based on smaller samples with weaker sig- 
nal/noise and lower dispersion data, only a minority of our 
sample (16%) appears to be due to genuinely young dust- 
reddened systems. We suggest the discrepancy may arise in 
part because of our choice of I — H as a colour discriminant 
rather than the coarser R — K. 

However, within our dominant population of established 
galaxies, we notice a great diversity in spectral properties. 
We classify these as pure E (27%), E+A (13.5%), E+e (23%) 
and 'mixed' E+A+e (13.5%) depending on the presence or 
otherwise of secondary features of star formation such as [O 
II] emission and post-burst Balmer absorption. We propose 
that these different subsets represent various manifestations 
of the same overall population seen at different stages as 
they accrete associated objects and via illustrative spectral 
modelling deduce that the bulk of the stars in the population 
formed well before z ~2. 

The abundance of the E+A and similar systems in our 
sample is an extremely interesting result. The numbers are 
quite high, in comparison with the local universe where they 
are very rare, and imply we may be seeing the fading light 
of a recent major star formation epoch. 

As earlier papers in this series have surmised, the abun- 
dance and luminosity distribution of this established com- 
ponent of EROs, when faded by passive evolution to the 
present epoch, cannot account for the population of local 
early-type galaxies. As many workers are deducing (Treu et 
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al 2005, Bundy et al 2005), continued growth and trans- 
formations are required over 0.5< z <1.5 to match the lo- 
cal distribution. However, as many authors have claimed, 
the evolutionary growth in the number of early-types since 
z < 4 is far less dramatic than that predicted in recent semi- 
analytical models. 

We finally refine earlier estimates of the spatial clus- 
tering of this population taking into account the improved 
redshift distribution in the inversion from the angular clus- 
tering seen in the overall LCIRS survey. We find a spatial 
correlation length, corrected to the present epoch assuming 
no evolution in comoving coordinates, of ro=6.5 ± 0.4 h~ 
Mpc, a signal comparable to but marginally less than that 
seen in present day ellipticals. 
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